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Analytical expressions and reduction formulas are developed for various indefinite integrals of the 
confluent hypergeometric functions. These integrals are of the type //(a, b, z)z , 'e az dz, where /is one 
of the two Kummer functions M(a, />, z) = \F\{a\ b; z) or E/(a, /;, z), with real or complex a. b, z and a, 
and real p. 
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1. Introduction 

Multiple integrals occur in many physical problems in the form where one of the integrations 
results in a higher mathematical function. Thus one is often faced with the task of evaluating in- 
tegrals involving higher functions. Whereas there exists a large body of literature for integrals with 
limits (0, oc) and (— », »), there are considerably fewer sources pertaining to indefinite integrals. This 
latter type is important, for example, in statistical mechanics for the integration of a distribution 
over a finite portion of velocity or momentum space and in electromagnetic theory for the integra- 
tion over a finite charge configuration. We note the appearance of only one major monograph with 
emphasis on indefinite and finite integrals of higher functions, that of Y. L. Luke [l] 1 on Integrals 
of Bessel Functions. This monograph covers a much wider range of material than its title suggests. 

In this note we investigate some indefinite integrals involving confluent hypergeometric func- 
tions of the type //(a, b, z)z p e az dz, where/is one of the two Kummer functions M (a, b, z) = \F \ (a; 
b; z) or U(a, 6, z) . Several recent applications of these integrals to studies on molecular electronic 
energies [2, 3] have been made. 

It is well known that the confluent hypergeometric functions contain many other higher func- 
tions as special cases. For the definitions and an extensive analysis of the properties of these 
functions, the reader is referred to the Handbook of Mathematical Functions [4] and the work by 
Slater [5]. 

We define the integrals to be considered in this note as follows: - 

M p (a; a, 6, z) = J M(a, b, z)z^e az dz, (la) 

U p (a; a, 6, z) = J U(a, b, z)z^e az dz, (lb) 

and the special cases 

® n (a, b 9 z) =M n (0; a, 6, z), (2a) 

V n (a, b,z) = U n (0; a 9 b,z). (2b) 



*An invited paper. This paper presents the results of one phase of research carried out at the Jet Propulsion Laboratory, California Institute of Technology, under 
Contract NAS7-100, sponsored by the National Aeronautics and Space Administration. Acknowledgment is made to Dr. M. Saffren for helpful discussions. 
**Present address: Jet Propulsion Laboratory, California Institute of Technology, Pasadena, Calif. 91103. 
1 Figures in brackets indicate the literature references at the end of tins paper. 
'Throughout this paper we omit the constant for all indefinite integrals. 
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Here a, 6, z, and a may be real or complex, p is real and n is a positive integer or 0. In section 2 
we discuss the function <& n and in section 3 the function ty n showing how both of these functions 
can be reduced to simple linear combinations of elementary functions and various hypergeometric 
functions. In section 4 we consider eq (1) with a = 0, and we conclude in section 5 with a brief 
discussion of the general integrals M p (a; a, b, z) and U p (a; a, b, z) of eq (1), relating them to some 
recent work by Babister [6]. In an appendix we specialize the parameters in M p (a; a, b, z) and 
U p (a; a, 6, z) to the Bessel function case, b= 2a, and, using the results of sections 4 and 5, display 
several interesting indefinite integrals for the Bessel functions. 

2. Analysis of <$> n (a, b, z) 

In this section we analyze the function <£«(#, 6, z) and show its reduction to elementary 
functions or to linear combinations of various hypergeometric functions. We use the following well 
known properties of the M(a, b, z) [4, 5] 

(l + a-b)M(a 9 6, z) - aM(a-+ 1, b, z) + (b-l)M(a, o-l,z)=0, (3) 

bM(a, 6, z)-6M(a-l, b 9 z)-zM(a 9 6+l,z)=0, (4) 

4~M(a, b 9 z)=rM(a+l 9 b+1, z); (5) 

az b 

and their equivalent integrated expressions (see footnote 2) 

(l + a-6)<Ma, 6,z)-a<I>»(a+l, M) + (6-l)* n (a, 6-1, z) = 0, (3a) 

60> n (a, 6, z)-6<Ma-l, 6, z)-4> M + ,(a, 6 + l,z) = 0, (4a) 

z n M{a, b 9 z)-/i$ n -i(a, 6, z) =T*n(a+l, 6+1, z), (5a) 

M(a 9 l,z)=7$ (fl+l,4+l,z). (5b) 

Equations (5a) and (5b) are obtained by integrating eq (5) by parts. We also note the useful 
expressions 

M(0, b 9 z) = l, (6) 

4> n (0, 6, z) =-£^j- (6a) 

Four different cases arise dependent on whether a is equal or unequal to an integer and b is 
equal or unequal to an integer. 

(i) a 7^ integer, b # integer 

Use of eq (5a) with a replaced by a — 1 and b replaced by b — 1 yields the reduction formula 

$»(aJ, Z )=j^z»M(fl^U-l,z)~nj^Mfl-U-l,^ (7) 

(a-1) (a-1) 

Repeated application of eq (7) generates ^>«(a, 6, z) as a linear combination of terms 
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z u ~ n 'M(a — n' — 1 , 6 — n f — 1, z), with ra' ranging from to n — 1, plus the lowest order term, 
<£>o(tt — /?, b — ju z), which is given by eq (5b). 

(ii) a # integer, 6 = integer > 

We first note that the integral, as well as M (a, 6, z), is undefined for b = — m. For 6 t^ 1 we can 
replace (b— 1 )4> M -i (a — 1 , 6 — 1, z) in eq (7) by its equivalent from eq (3a) with a replaced by a—\ 
and n replaced by // — 1 . , to vield a reduction formula analogous to eq (7) 

<£> B (q, b, z)= (6 ~ 1 j z w M(a-l,6-l,z)-n4> w -i(a, 6, z) +/i j fl ~^ <D w -i(a- 1, 6, z). (8) 

Repeated application of eq (8) generates for <J> n a linear combination of M's plus the lowest order 
terms Oo(a', b' , z) which we obtain from eq (5b). For the case b= 1 we see that eq (8) is still valid 
if we recognize from eqs (3) and (4) that 

Lim (b-l)M(a-l, b- 1, z) = (a- l)zM(a, 2,z). (9) 

b-* l 

Substituting eq (9) into eq (8), we have the reduction formula 

<£ n (a, 1, z) = z" + 1 M(a, 2, z) -7i[4> n -i(a, 1, z)-* n -i(a-l, 1, z)], (10) 

with the lowest order term being given by 

4>„(a, l,z)=zM(a, 2, z). (11) 

One can easily show by mathematical induction that the general result for the case b= 1 is 



<D„(a, 1,2) =2-' Y (-l)"(' l )7-XTT M («' m + 2 ^) 
£. \m) (m+1) 



(12) 



which holds equally well when a is an integer. 

(iii) a = integer, b # integer 

For a > n + 1, eq (7) can be used to reduce ^(a, 6, z) to a series of M's plus the lowest order 
term <f>o(a — n, b — n,z) given by eq (5b). For a=n+l and a = ft > 0, repeated application of eq (7) 
again reduces <Pn to a series of M's plus the lowest order terms 4>o(l? b — n,z) and $i(l, 6 — n-h 1, z) 
respectively. These are related through eq (4a) with a= 1 and ra = and eq (6a) by 

*i(l, 6-Ti + l, z) = (b-n) [<D (1, 6-n, z)-z]. (13) 

The primitive function ^oQ, b' , z), not reducible any further, can be written [1, p. 5] 

0>o(l, 6',z)=z 2 F 2 (l, l;2,6';z). (14) 

For the special case where b' is an integer, the 2F2 function can be reduced to incomplete gamma 
functions and exponential integrals as will be seen in case (iv). The condition n > a > is easily 
handled by repeated use of eq (4a) with n replaced by n — 1 and b replaced by 6 — 1 

$„(a, 6,z) = (6-l) [<&,,_, (a, 6-l,z)-*„-i(a-l, 6-1, z)], (15) 

until the reduction generates either of the lowest order terms <lV(0, 6, z) given by eq (6a), 4>o(«' ^ 1» 
6', z) of eq (5b), or <t>o(l, 6', z), the primitive given in eq(14). Lastly, for the case a = — m, m=0, 1, 
2, . . ., since M(— m, 6, z) is a polynomial of m+ 1 terms we can either employ eq (5a) backwards 
until we arrive at <t> n + m(0, b + m, z) , given by eq (6a), or directly integrate to obtain 
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**(-**>!>>*)= 2 (-W {' k ) (b) ; {n+k+1 y (16) 

where 

(b) k =r(b + k)ir(b). 

(iv) a = integer, b = integer > 0. 
With a 5? b, repeated use of eq (4) generates 

M(b + m,b, *)-*£$)■£;, (17) 

since 

M(b,b,z)=e z , (18) 

which holds even when b t^ integer. Integration of eq (17) obviously yields 

*.(A + M . *.*>-* (T) cfe / 



(19) 



k=0 

For the case b > a > 0, we have from eq (3a) with b = a + m, the recursion relation 

(1 — m)<$> n (a, a + m, z) = «<£„(« + 1, fl + m,z)-(a + wi-l)$„(a, a + m—l, z) . (20) 

Equation (20) is repeatedly applied until the right hand side contains only terms of the form b = a+ 1. 
Using eq (4a) with b replaced by a and n replaced by n — 1 , we have the reduction formula 



O n (a, a+1, z) =a I z n ~ x e z dz — a<& n -i(a— 1, a, z). 



(21) 



The reduction by eq (21) is continued, generating a series of terms fz k e z dz plus the lowest order 
term. If a ^ n, the lowest order term is O„-„(0, 1, z), which is given by eq (6a). With a = n+l, 
the lowest order term is 

*o(l,2,*)=J"(e»-l)^, (22) 

or equivalently 

* (1, 2, x)=Ei(x)-\nx x>0, (23a) 

<J>o(l, 2,-z)=-E x (z)~ Inz |argz| < *r, (23b) 

where £\ and £7 are the exponential integrals [4, Chap. 5]. Finally, if a > n+ 1, the lowest order 
term is <Po(a — n, a — nH- 1, z). Using eq (5b), we can express this in terms of M(a — n— 1, a — n, z) . 
Since 

M(n, n + l,-*)=^ [1-c'ViW], (24) 

where e w (#) = V tj; we have 
fc=o *• 

g>o(a-/i,a~/i+l,z) = (-l)«-^ (a ~^l a ~r~ 2)! [l-e g e ff - w - 2 (-z)]. (25) 
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We note that since 




yU, x) =(i l x a e-*M(l, 1 + a, x) = a l x"M(a, l + a,-x), 


(26) 


where y(a, x) is the incomplete gamma function [7], the 4>o(l, b\ z) of eq (14) can al 


teniatively 


be written as 




<t>o( 1, 6\ z) = (6'- 1) f y(6' - 1, zjz'-^fife, 


(27) 



whereas <J>oU/, 1 + a, — x) can be expressed as 



For g/1, this reduces to 



<!>,,( a, 1 + a, —x) = a I y(a, x)x a dx, 



O (a, a+1, -x) = y^— [^^^(a, %) +e"- r ]. 



(28a) 



(28b) 



For a= 1, eq. (28a) reduces to eq. (23b), and with a= integer > 1, eq (28b) reduces to eq (25). 

The last case to be considered in this section is b > a where b > and a=—m. This case is 
immediately obtainable from eq. (16) since that equation obviously holds for b= integer > 0. 

3. Analysis of ^,,{(1, />, z) 

Analogous to the results of section 2, the function *P n (a, b, z) can be reduced to elementary 
functions or to linear combinations of various hypergeometric functions. Here, too, four different 
cases, a equal or unequal to an integer and 6 equal or unequal to an integer, must be considered 
for the complete reduction of M',,. We commence the analysis with the definition of U(a, b, z) [5] 



U(a, 6, z) 



sin 77 b 



M(a, 6, z) 



h M(l + a-b, 2-6, z) 



l(l + a-b)l(b) r(a)r(2-6) 

and the well known properties [4, 5] 

U(a, 6, z)-aU(a 4-1, 6,z)-tf(a, 6-1, 2) =0, 
(6-a)J/(a, 6,z) + f/(a-l, 6, z)-zU(a, 6 + l,z)=0, 



<b 



f/(a, 6, z)=-atf(a+l, 6 + 1, z). 



We also note the integrated form of eqs (30) — (32)- 

*„(a, 6, z) -a¥„(a+l, 6, z)-"* n (a, 6-1, z) =0, 

(6-a)¥ n (a, 6, z)+^„(a-l, 6, z) -* w+ i(a, 6 + l,z)=0, 

z"f/(a, 6, z)-/i* n -i(a, 6, z) =-aV n (a+ 1, 6+1, z), 

f/(a, 6, z)=-a^ (a+l, 6+1, z). 



(29) 

(30) 
(31) 

(32) 

(30a) 
(31a) 
(32a) 
(32b) 



Equation (32b) is simply obtained from eq (32a) by setting n = 0. Other properties of the U(a, b, z) 
used in this analysis are 

f/(0, 6, z) = l, (33) 



M'„(0, 6,z)=z w+1 /(* + l), 

c/(a, a+1, z)=z~", 

£/(a, 1 — n, z)=z"c/(a + rc, 1 + n, z), 
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(33a) 
(34) 
(35) 



the last equation being a special case of Kummer's transformation. One final property of the 
U{a, 6, z) , as opposed to the M(a, b, z) , is that they are defined for all values of a and b. 

(i) a 7^ integer, b 7^ integer 

Using eq (32a) with a replaced by a — 1 and b replaced by 6 — 1, we have 

^„(a, b, z)= * r [z n U(a-l 9 b- 1, z) -n9 n -i(a- 1, 6- 1, '*)]. (36) 

v t ^/ ♦ 

Repeated use of eq (36) reduces ^ n to a series of terms z n ~ n 'U(a — n' — \,b — n' — \,z) plus the 
lowest order term, which by eq (32b) is given by 

¥o(a\ b\ z)= l ,. U(a'-l, b'-l, z), (37) 

(1-a ) 

or equivalently, from eq (30) with a replaced by a — 1 

*o(a\ b\ z) = l ,. U(a'-l, b\ z) + U(a\ b\ z). (38) 

(1-a ) 

(ii) a # integer, b — integer 

For b ^ 1, if we replace both terms on the right hand side of eq (36) by their equivalents from 
eq (30) and (30a) with a replaced by a — 1, then we have the reduction formula for NPn 

¥„(a, b, z) = [z n {U(a-l, b, z) + (l-a)U(a, fc, z)} 

-7i{¥ tt -i(a-l, fc,z) + (l-a)^„-i(a, 6,*)}]; (39) 

the lowest order term is given by eq (38). If b ^ 0, we use the integrated form of Kummer's trans- 
formation, eq (35), to yield 

*„(a,-77i,z)=*„+ m +i(a + 77i+l, tti + 2, z), (40) 

which can obviously be handled by means of eqs (39) and (38) above. 

(hi) a = integer, b t 4 integer 

If a = — tfi, 77i = 0, 1, 2, . . ., then from eq (31a) with a = — ra-f 1, we have the recursion 
relation 

Vni-m, b, z) = ¥» + i(- TTi+l, 6+1, z)-(6 + 77i-l)¥ n (-/7i + l, 6, z), (41) 

which reduces to ^'(0, &', z) of eq (33a) upon repeated usage. An alternative approach is to note 
that since 

U(-m, b,z)={-l) m m\W-^{z), (42) 

where Lfa)(z) is the associated Laguerre function [4, chap. 22], direct integration yields 

•^»..)-(-.)^ I l(-:!; 1 )^^ iy . 
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For a > 0, three distinct situations arise. If a >/i + l, use of eq (36) reduces ^(a, b, z) to 
a series of f/'s plus the lowest order term v J / (a — n, b — n, z) given by eq (37) as in case (i) above. 
With a = n or a =/i-f L, repeated application of eq (36) reduces ^ n to a series of t/'s plus the low- 
est order terms ^o(l, 6 — n, z) and ^i(l, 6— /* + 1, z) respectively. These are related through 
eq (31a) with n = and a= 1 by 

¥i(l, 6' + l,z)=z+(6'-l)* (l, 6',*). (44) 

Furthermore, from the definition of U{a, 6,z), eq (29), with a = 1, we see that 

•¥ Q,6',z) = r(6'-l) \j£rdz- l 3>,(l,6',z), (45) 



where the first term in eq (45) is an incomplete gamma function and the second term is precisely 
the primitive function discussed in section 2, eq (14). We remark that for the case b = integer, 
eq (45) reduces to powers, logarithms and the exponential integral as will be seen in (iv) below. 
Finally, for a < n, use of eq (31a) with n replaced by n—\ and b replaced by b — 1 leads to the 
reduction formula 

V n (a, &,z) = (6-a-l)¥„-i(a, 6-1, z) + ¥ tt -i(a-l, 6-1, z). (46) 

Continued application of eq (46) results in the lowest order terms ^,'(0, b' , z) given by eq (33a) 
and tyn'in' , b\ z) , the case a = n discussed above. 

(iv) a = integer, b = integer 



Corresponding to eq (26) for the representation of the incomplete gamma function, y(a, x), 
in terms of M , we have 

T(a, x)=x a e-*U(l, l + a,x) = e-*U(l-a, l-«,x). (47) 

Therefore for the case a =b = — m, m = 0, 1 , 2, . . ., we have 

^„(-m,-w,z)= I z n e*r(m+l, z)dz. (48) 

Replacing r(m-\-l, z) by its expansion as a truncated exponential and integrating termwise, 
we obtain 

m ] 7 n+k+\ 

*.(-„,-„,,)-„! g - ? - f _ (49) 

For a = b= z -\-m, we have a two-step reduction procedure. First, from eq (30a) with a= m and 
6= Hi-f 1 we have 



1 f 1 

*n m+1, /n+1, z)=— z"- m dz *«(m, in, z) , 

m J m 



(50) 



where we have used the integrated form of eq (34) to eliminate ty n (m, m+ 1, 2). This reduction 
is continued until ^„ ( 1 , 1 , z) is generated. Second, noting that 

1/(1, l,z) = e*£i(z), (51) 
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multiplying by z" and integrating by parts, we obtain the reduction formula 



¥„(1. l,z)=z n e*E 1 (z)+ ± — n¥„-i(l, 1, z) 



(52) 



Equation (52) is applied n times until ^o(l, 1, z) is generated which, by integration of eq (51), is 
obviously given by 



*o(l, l,z) = lnz+e*£i(z). 



(53) 



This same result is obtained from eq (47) with a chosen as 1 — m and the replacement under the 
integral sign of F(l — /ra, x) by x x ~ "' E,„( x) . 

For b > a, if a ^ 0, eq (41) is repeatedly applied until a series of terms ^IV(0, b\ z) is obtained 
which is given by eq (33a). With a > 0, specifically choosing a = ra+land & = m + ra'-h2 in eq 
(30a), we have the recursion relation 



¥„(m+l, ira + m' + 2, z) = (m+ l)^„(ra + 2, ni + m' + 2, z) +^»(m+l, ro + m' + l, z). 



(54) 



This recursion is repeated until the right-hand side contains only terms with 6 = a + l, given by 
the integral of eq (34). The general result is then simply 



¥h(iw+1, m + m' + 2,z) : 



1 



(m + A) 



'/*"■ 



l dz. 



(55) 



Finally, the case a > b is either treated by the reduction formula eq (41) for a ^ until ^n'(0, 
b\ z) of eq (33) is obtained, or for a > 0, by repeated use of eq (31a) until a = or a = b=-\- m 
(considered in eqs (50)— (53) above) is reached. Often one can utilize eq (35) to reduce the parameters 
to lower orders before starting the recursions. As a simple example of the case a> b, we present 
the result for * 2 (4, 2, z) 



¥,(4, 2,z)=%(3, 0,z) = j2 [(l + *)-*(2 + z)e*£i(z)]. 
4. M p (0; a, 6, z) and E/ p (0; a, 6, z) 



(56) 



Consider first the integral M,,(0; a, 6, z) . When the series for M(a, b, z) is integrated termwise, 
we find that 



*,<*..,.„- J $jXIfe±gil, -£.«,+!..,,+..»;.>. 



(57) 



provided that 6^0, —1, —2, . . ., and p# — 1, — 2, . . . . We note that for p = — m, with 

m=l, 2, . . ., 



M_ m (0; a, 6, z) = ( ^; m * 3 F 3 (tt-fm, 1, 1; 6+ m, 2, l + m;z) 



(b) m m V 



m y 2 (a)t 1 



1 



(6)**!(i»-i-l) 



' + 



(a). 



] 



(6) m -i(m-l) ! 



In z. (58) 



From the definition for £/(a, b, z) in terms of M(a, 6, z), eq (29), we have 

1 



U p (0; a, 6, z)=-t 



sin 7r/> 



r(l-r-a-fe)r(6) 



M p (0; a, 6, z) 

1 

r(a)r(2-6) 
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M„ +1 _ 6 (0, l + a-6, 2-6, z) 



(59) 



Contiguous relations for the generalized hypergeometric functions, ,F <h are given by Rainville 
[8]. From the special cases for r = q = 2, we obtain 

aM p (0; a + 1, 6, z)=z*> +1 M(a, 6, z) - (p + l-a)Af,(0; a, 6, z), (60) 

(6-l)Af p (0; a, 6-1, z) =z" +1 M(a, 6, z) - (p + 2-6)M p (0; a, 6, z), (61) 

and for the U p (0; a, 6, z), from eqs (59), (60), and (61), we have 

a(l + a-b)U p (0; a -hi, 6, z)=z" M [/(a, 6, z) - (p + l-a)f/ y ,(0; a, 6, z), (62) 

(6-a-l)f/ p (0; a, 6-1, z) =z" + l f/(a, 6, z) - (p + 2-6)t/ p (0; a, 6, z). (63) 

Finally, by integrating eq (1) by parts, we have the recursion relations 

aM p+l (0; a+1, 6 + 1, z) =b[z» +l M(a, 6, z) - (p + 1)M P (0; a, 6, z)], (64) 

af/ p+ i(0; a-f 1, 6 + 1, z) = (p + l)£/ p (0; a, 6, z)-z» +l U{a, b, z). (65) 

For p = integer n, n = 0, 1, 2, . . ., M„(0; a, 6,z) = 0„(a, 6, z) and £/„(0; a, b,z) = ^„(a, 6,z), 
the cases discussed in sections 2 and 3 respectively. For/; # integer, four cases arise for which the 
integrals can be evaluated in terms of elementary functions or in terms of AT s or f/'s respectively. 
The first case is the trivial one with a = 0. We obviously have, by eqs (6) and (33) 



M p (0; 0, 6, z) = U P (Q; 0, 6, z) =-^. (66) 

/> + 1 

The second case, a = p + 2 + m, m = 0, 1, 2 is reduced by repeated use of eqs (60) and (62) 

to the lowest order terms a = p + 2, where 

M p (0; p + 2, 6,z)=^|M(p+l, 6,z), (67a) 

U p (0; p + 2, b, z) = — f '' tf(p + I, 6, z), (67b) 

(p+ L)(p + 2 — 6) 

which are directly obtainable from eqs (60) and (62) by setting a = p-\- 1. Thirdly, if 6 = p+l — m, 
then use of eqs (61) and (63) leads to the lowest order terms 6 = p+ 1 , where 

M p (0; a,p+l, z) =-^-rM(a,p + 2, z), (68a) 

U p (0; a,p + l, z)=(-~^^W(a,p + 2,z), (68b) 

which we find from eqs (61) and (63) by setting b = p-\-2. Lastly, for b — a = ± integer, use of eqs 
(60) and (61) reduces M p to the case a= b for which 

M p (0; a, a,z)=-^j-M(p+l,p + 2,z), (69a) 

whereas, if 6 — a=l, 2, . . ., repeated application of eqs (62) and (63) reduces U p to the case 
b = a + 1 for which 



U p (0; a, a+l,z)= I z»-"<Iz. 



(6%) 



We note that if none of the above conditions are satisfied, then use of eqs (60)-(65) will reduce the 
integrals to a series of similar integrals with p, a, 6 in the interval (0, 1). Further analytical reductions 
are not available. 
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5. The General Integrals M p (a; a, ft, z) and U p (a; a, 6, z) 
Let us consider the differential equation 

z, B +{b ~ z) fr a ^ f{z) - (70) 

For/(z) =0, we have the well-known confluent hypergeometric equation for which M(a, 6, z) and 
£/(a, 6, z) are solutions. Babister [6] defines two nonhomogeneous confluent hypergeometric 
functions 6<r(a, 6; z) and A p ,cr(a, 6; z) which are particular solutions of eq (70) for (i)f(z) =z iT ~ l 
and (ii) f{z) =exp (z)z <T ~ 1 respectively. He further discusses series expansions, contiguous rela- 
tions, recurrence formulas, differential properties and several definite and indefinite integrals 
for 6 and A. 

We shall presently express M p and U p in terms of the functions investigated by Babister. 
Consider the following special cases of eq (70), 



and 



z^+(b-z) ^-ay=exp (pz)z-- 1 , (70a) 

dz l dz 

z^+(b-z)^-au=0, (70b) 

dz- dz 



with particular solutions y= A p ,<r(a, b; z) and u = M(a, 6, z), respectively. Multiplying eq (70a) 
by ue~ z z b ~ x and eq (70b) by ye~ z z h ~\ and subtracting, we obtain 



A 

dz 



\-*\M{a, 6, z)— A p ,<r(a, 6; z)-A p , (T {a, 6; z) — M(a, 6, z)j 

= e <P-i)z 2 !>«T-2M(a, 6, z). (71) 
Integrating both sides of eq (71) and setting b + <r — 2 = p and p = a+l, we have 

M p (a; a, 6, z)=z b e~ z M(a, 6, z) — A a +i, p -b+2(a, 6; z) 

-A a+ i, p _ 6+2 (a, b;z)^- z M(a, b, z)\ (72) 

From the differential properties of M(a, 6, z) , eq (5), and A p , (r (a, 6; z) [6, p. 135] 

-^A p , C7 (a, 6;z)=(<r- l)A Pf<r -i(a+l, 6 + 1; z) +pA p , <,(<*+ 1, 6 + 1; z), (73) 

dz 

we obtain an alternative expression for M p (a; a, 6, z) of eq (72) 
M p (a; a, 6, z)=z /> e- 2 ' 



(a+l)M(a, 6,z)A«+i, p- 6 + 2 (a +1, 6+1; z) 
+ (p — 6+l)M(a, 6, z)A a +i,p-6+2(a+l, 6+1; z) 



■| M(a+1, 6+1, z)A a +i,p-6+2(fl, 6; z) 



(74) 



Contiguous relations for the A are given by Babister [6, p. 135]. Two interesting special cases of 
eq (74) are obtainable, when a = and when a = — 1 due to the relations [6, p. 134] 

A ,o-(«, 6; z) = Said, 6; z), (75a) 

Ai, ff (a, 6; z) = e- i ° r < T e z (r (b-a, 6, e /7r z). (75b) 
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The first is obtained from eq (74) with a = and eq (75b) and gives an expression for the A/ p (0; a, 
6, z) of section 4 in terms of 0's and M's 



M ; ,(0; a, 6, z) = z*e ,,r(ft -^ 



M(a, 6, z)dp- b +2(b-a, 6 + 1; e /7r z) 
(p-6+l)M(a, 6,z)(?p- 6+ i(6-a, 6+1; e /7r z) 



-- M(a+1, 6+1, z)6 p - b+2 {b-a, 6, e /77 z) 



Furthermore, from eq (57) and from the following expression [6, p. 121] 



daia, b; z) 



o-(o-+6-r 



2^2(1, cr+a; cr+ 1, cr+6; z), 



we obtain 

2^2(7^+1, a; p + 2, 6; z)=M(a, 6, z) 2/2U, p-a+1; p-6 + 2, p + 2; -z) 

+ £ ^4t 7 . =r M(a, 6,z) 2 F 2 (l,p-a + 2; p-6 + 3,p + 3; -z) 



p + 2 (p-6 + 2) 



(76) 



(77) 



b (p-6 + 2) 



— M(a+1, 6+ I, z) 2 F 2 (1, p-a + 2; p-6 + 3, p + 2; -z). (78) 



The second case of interest is obtained from eq (74) with a = — 1 and eq (75a), and 
is given by [6, p. 131] 



M p (-1; a, 6,z)=zV-* 



(p-6 + l)M(a, 6, z)0 p -&+i(a + I, 6+ 1; z) 

-|M(a+l,6+l,z)fl p _ 6+2 (a, b; z)\ (79) 



M p (— 1; a, 6, z) can be related to the M ;> (0; «, 6, z) discussed above and in section 4, if we use 
Kummer's relation [4, 5] 



e~ z M(a, 6, z)=M(b — a, 6, c' v z); 



thus 



M p (-l; a, 6, z) = e- /7r ^ +1) A/„(0; 6 -a, 6, e^z). 
The condition 6 = a in eq (81) gives 



(80) 



(81) 



M p (— 1; a, a, z)- 



zp+ 



p+V 



(82a) 



corresponding to the first case considered in section 4, eq (66). The integral of eq (80) in addition 
to eq (5b) gives 



M (-l; a, b, z) 



b-\ 
a+l-b 
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e- z M(a, 6-l,z). 



(82b) 



Cases two and three of section 4, eqs (67a) and (68a), allow for the analytical integration of M p (— 1; 
a, 6, z) when b = a + p + 2 + m and 6=p+l— m for m = 0, 1, 2, . . ., and reduce respectively 
to the lowest order terms 



T J)+] 



M p (-1; a, a + p + 2, z)= J =— e~*Af(a+l, a+p + 2, z), 



Af ;> (-1; a,p + l,z) = 



7 p+i 



P + l 



e"*M(a+l,p + 2, z). 



(82c) 
(82d) 



From case four of section 4 we can obtain the integrals M p (—l; ±n, 6, z) for n = 0, 1, 2, 
Specifically, for ra positive we have 



M p (-1; n, 6, z)=zi J + l e~ 



^ b rrlk z^Mihp + ^z) 



(p-6 + 2)„p + l 



+ X<- 1 )*7^: 



(b-n) f 



(p-b + n+l-k) k + 



M{n-k, b, z) 



(82e) 



For /i negative, rather than using the reduction formula eq (60), it is simpler to integrate term by 
term the finite polynomial M( — n, 6, z) . 

The analysis proceeds in a straight forward manner for the U p (a; a, 6, z) if we use the definition 
of U(a, 6, z) in terms of M from eq (29). We then have 



U t ,(a; a, 6, z) '■ 



sin Trb 



1 



r(l + a-i)r(A) 



Mp{a\ a, 6, z) 



1 
Y{a)Y{2-b) 



M t ,- b + i(a; a — 6+1, 2-6, z) 



We first note that for b = a + 1 



U p (a; a, a + 1, z) = M p - a (a; 0, 1-a, z) = \zP~ a e a ' 



dz. 



(83) 



(84a) 



We also obtain from eq (83) for a=— 1 and eqs (26) and (82e) the case b= a 



U p (- 



1 



1; a, a, z) =-— [z" + T(l -a, z) + y(p- a + 2, z)]. 



(84b) 



Analogous to eq (82b) we find that by repeated use of eq (31), U„(— 1; a', 6', z) is reduced to a 
series of terms £/<)(— 1; a, 6, z), given by 

J7o(— 1; a, 6, z) =-e-*U(a, 6-1, z). (84c) 

Corresponding to the cases 6 = a + p + 2 + m and 6 =p + 1 — m we have 

m m i 

t/p(-l; a, a+p + 2 + m, z)=-z^ + 1 e-^ V (-1)* — — U(a+l + k, a+p + 2+/n, z), (84d) 

^o (m — /r)! 



£/ P (-l; a, p+l-m, z) =-^ +1 e" 2 ]£ 



e () (^-A)! 



f/(a + l + A, p + 2-m + A:, z). 



(84e) 



We note that for b = ±n, care must be taken in the use of eqs (59) and (83). In this instance, 
the limiting forms for these equations must be used. For a detailed discussion of this point, the 
reader is referred to Erdelyi et al. [9]. 
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6. Appendix: Reduction to the Bessel Function Case 

If 6= 2a, then flf(a, b, z) and U(a, />, z) are related to the modified Bessel functions / and K 

M\v+~, 2i>+1, 2z) = r(i^+l)2^- , 'e af /„U), (Al) 

tf(*+jL 2*>+l,2z) = tt- , / 2 (2z)- , 'c^U). (A2) 

Specializi. g the parameters in eq (57) to the case b = 2a and substituting from eq (Al), we have 
the indefinite integral 



I z»e ± *I v (z)dz = 1 , XT^Tv XTT 2 " 1 ^ U + p+1, i>+^;/x + i/ + 2, 2i>+l;±2zV 



(A3) 



in agreement with Luke [1, p. 96). Equations (67a) and (82c) give special cases of eq (A3) with 
/ll=— 3/2 and p= s p+ (3/2) whereas eqs (68a) and (82d) give special cases of eq (A3) with p,=v=p/2. 
Equations (67b) and (84d) of the text give rise to 



J z- 3/ VK P+ (3/2)Wtfe=- V2 



(2z)>" 
(p + l)(p + 2) 



^ [(/.-f2-z)/C / , + ,, /2) (z)-z^ +( , /2) (z)], (A4) 



(p + l)(p + 2) 
jz-%^ + (3/2)(# = - V2^ (2z)"+V 2 *tf(p + 3, 2p + 4, 2z) 



z- ] ' 2 e 



— [(p + 2+z)K pH zj 2 )(z) -zK p+ w2)(z)]; (AS) 



(p + l)(p + 2) 

whereas eqs (68b) and (84e) yield 

\ z v e*K v (z)dz = y/n % Z ;* *z v U(v+U2, 2i/ + 2, 2z) 
J (21/4-1) 



(2i> + 1 ) 

J z I 'e- 2 J K,(z)dz = -V^(2z)^ t ' +1 e- 2 ^(^ + 3/2, 2^ + 2, 2z) 

=^^ ) [K v (z)-K v + l (z)], (A7) 

in agreement with the Handbook of Mathematical Functions. 3 We note further that eq (84b) gives 

J e-*z-*K v (z)dz=j^zrne-*z i ->'[K v . 1 (z)-K v (z)] (v*l/2). (A8) 

From eqs (74) and (78) we obtain indefinite integrals of the product of powers of z, Bessel 
functions, and exp (<xz) expressed in terms of the A's and 0's of section 5. Equation (79) leads to 



3 See eqs. 11.3.15) and (11.3.16) in 1 1, chap. 11]. Luke [1, p. 97] gives an equivalent form for (A6) and (A7)when v = (). 
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an interesting identity involving 2/Vs 



(fWlT^+l)]- 1 2 F,(^+l/2, n + v+1; 2v+\, jjl + p + 2; -2z) 

= Ir(z) 2 F 2 (l, /ll+|; M-^+l, fji + v+2; 2z\ 



(p-v+i) L 



/rU)+/ F+ l(z) 



2 F 2 (1, /x+|;/Lt-v + 2, m + ^ + 2;2z). (A9) 



Finally, from eq (82b), we obtain a decomposition formula for 2 F 2 (a, b; c, of; z) when a=l, c = 2, 
and d=2b. We then have 4 

z 2 F 2 (l, a + 1; 2, 2a + 2;z)=- (a+l/2)[M(a, 2a + l,z) — 1], (AlOa) 

for a # — A/2, A; = 0, 1, 2, . . ., or in terms of Bessel functions 

2z 2 F 2 (l, p + 1/2; 2, 2*+l; 2z) - ^ [ (2/z)*ze*I> + l){J„_i (2) -/„(*)} -2v]. (AlOb) 
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4 Equation (A 10a) is a particular example of the more general result derivable from the series representation oi p F q . 
WVmU, a, + 1; 2, *,+ !; *)- [„f,(o„; b,; j)-l] (f\ a,) I I f[ b,\ 
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